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We use a microscopic density functional theory based on Wertheim's first order thermodynamic 
perturbation theory to study wetting behavior of athermal mixtures of colloids and excluded-volume 
polymers. In opposition to the wetting behavior of the Asakura-Oosawa-Vrij model we find the 
polymer-rich phase to wet a hard wall. The wetting transition is of the first order and is accompanied 
by the prewetting transition. We do not find any hints for the layering transitions in the partial 
wetting regime. Our results resemble the wetting behavior in athermal polymer solutions. We point 
out that an accurate, monomer-resolved theory for colloid-polymer mixtures should incorporate the 
correct scaling behavior in the dilute polymer regime and an accurate description of the reference 
system. 
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I. INTRODUCTION 

Mixtures of nonadsorbing polymers and colloidal par- 
ticles often exhibit rich phase behaviou^. For certain 
polymer-to-colloid size ratios q = R g /Rc (where R g is 
the radius of gyration of the polymer and Rc is the radius 
of colloid) entropy-driven effective interactions may lead 
to stable colloidal gas, liquid and solid phases even if all 
bare interactions are purely repulsive 3 -. A simple theoret- 
ical model giving an insight into this phenomenon is the 
Asakura-Oosawa-Vrij (AOV) model of colloid-polymer 
mixtures in which the polymers (modeled as spheres) are 
ideal and can overlap freely, but the polymer-colloid and 
colloid-colloid interactions are of the hard sphere type^. 
Effective attractive interactions in such systems arise due 
to a tendency to decrease the volume excluded to the 
polymer coils£. 

The AOV model has attracted much attention due to 
its simplicity^, however when comparing to experiments 
the agreement is only qualitative. This is due to the 
fact that the real polymers are nonideal. The incorpo- 
ration of the polymer nonideality can be tackled at the 
monomer-resolved&2ii2iii or coarse-grainedi2iiii4 levels 
of description. The latter technique is particularly useful 
in the so-called "colloid limit" (q < 1), in which the full 
mixture is well described by invoking the effective, pair- 
wise depletion potentials 15 . However, when the polymer 
dimensions are larger than the size of the colloid particles, 
the resulting two-body effective potentials may be insuf- 
ficient to correctly describe the underlying mixture and 
the incorporation of the higher-order, many-body terms 
becomes necessaryiSiii. 

In addition to many industrial applications colloid- 
polymer mixtures offer a convenient tool to study im- 
portant fundamental concepts that are often difficult to 
investigate in simple fluids, such as real space observation 
of the thermal capillary wavesi^, the capillary length^, 
the interfacial width^S and the interfacial tension of the 
fluid-fluid interfaco 21 i 22 i 23 . Wijting et alM^ studied 
the behavior of colloid-polymer mixtures close to a wall. 



They found a wetting transition to a state in which the 
colloid-rich phase wets completely a nonadsorbing pla- 
nar hard wall. These findings qualitatively agree with 
the results of density functional theory 2 ^ and computer 
simulations 2 ^ carried out for the AOV model. Aarts et 
al^— used the Cahn-Fisher-Nakanishi approach to study 
the wetting transition in mixtures of colloids and ex- 
cluded volume polymers and pointed out some subtleties 
associated with the precise measurements of the contact 
angled. 

Recently Paricaud et al&L examined the bulk phase 
behavior of mixtures of hard-sphere colloids and 
excluded-volume polymers within the framework of 
Wertheim's first-order thermodynamic perturbation the- 
ory (TPTI)2i. They found a demixing transition into 
the colloid-rich and polymer-rich phases provided that 
the ratio of the colloid diameter to polymer segment di- 
ameter is large enough. The TPT1 description employs 
the microscopic description of the polymer, consequently 
it is treated on equal footing with the colloid particles. 
This theory provides thus a straightforward procedure of 
taking the attractive/repulsive interactions into account. 
However, the Wertheim TPT1 approach is not free of 
its own deficiencies. We note here that incorporation of 
the solid-fluid equilibrium into the framework, while in 
principle possible 3 ^, is difficult to execute for the colloid- 
polymer mixtures. Moreover, Boublik alM^ pointed 
out that within the TPT1 theory the second virial co- 
efficient scales quadratically with the chain length M in- 
stead of oc M 3l/ , although it is possible to improve this 
deficiency 3 ^. Finally, the value of the Flory exponent, 
v, resulting from the TPT1 approach is 0.5 rather than 
0.588 30 . Despite these shortcomings, the TPTI approach 
to colloid-polymer mixtures provides a welcome depar- 
ture from the well-studied AOV model. 

In this work we focus on the surface phase behavior 
of model colloid-polymer mixtures. We employ nonlo- 
cal density functional theory that was used previously 
to investigate liquid-liquid interfaces of athermal mix- 
tures of colloids and excluded-volume polymers 3 ^. The 
density functional framework, proposed originally by Yu 
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FIG. 1: Bulk phase diagram for a colloid polymer mixture 
with d = 10 and M — 100. The main figure shows the binodal 
(solid line), the spinodal (dashed line), and the bulk critical 
point (open circle) evaluated in the polymer segment mole 
fraction-total density plane. The open circle denotes the bulk 
critical point. The inset shows the binodals resulting from the 
TPT1 (solid line) and free volume (dash-dotted line) theories 
evaluated in the colloid-polymer packing fraction representa- 
tion. 



and Wu2£, uses the TPT1 ideas and yields in the limit of 
bulk systems an equation of state identical to that from 
RefiSI. Since the DFT framework has already been de- 
scribed in detail2Si2i we recall only the basic points of 
this approach. 



II. THEORY 

We model the colloids as hard spheres of diameter 
ac and the polymers as chains composed from M tan- 
gentially bonded hard-sphere segments of diameter aps- 
The hard-sphere monomers forming the chains are freely 
jointed i.e. they can adopt any configuration as long as it 
is free of the intermolecular and intramolecular overlap. 
Within the DFT approach the grand potential of the sys- 
tem, f2 is a functional of the local densities of polymers, 
pp(R) and colloids, pc{ r ) 

n[p P (R),pc(r)} = F[p P (K),pc(r)}+ 

JdRpp(B.)(v£>(1l)- - ' 



■ Vp) 



dvp c {v){V^{v) 



Vc) 



(1) 



In the above V}^ t \r), p c 



and pp are the ex- 



ternal and the chemical potentials for colloids and poly- 
mers, respectively. R = (ri, r 2 , • • • , ym) denotes a set 
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FIG. 2: The excess polymer segment adsorption isotherm 
Yps for a colloid-polymer mixture with d = 10 and M = 100, 
calculated along path I from Fig. (upper panel) and along 
path II from Fig. Q (lower panel). The dotted lines denote 
bulk coexistence. 



of monomer coordinates. The free energy of the sys- 
tem F is a sum of the ideal and excess contributions, 
F = Fid + F ex . The ideal part of the free energy is known 
exactly, (3F ld [p P (R), Pc (r)] = JdRp P (R)V b {R) + 
JdRp P (R)[\n(pp(R)) - 1] + Jdrpc(r)[ln(pc(r)) - 1]. 
The total bonding potential, Vb(R), is represented 
as a sum of the bonding potentials vj, between the 
monomers, Vb(R) — Yli=i 1 v b{\ r i+i — r i\), and satisfies 
exp[-(3V b (R)} cx n^I^d^+i - r *l - °ps)- Further- 
more, it is assumed that the excess free energy is a func- 
tional of the local density of colloids and average seg- 
ment local density defined as pps(r) 
ZfijdRSir-r^ppiR). 



J2tiiPPsA r ) 



Within the approach of Yu and Wu2£ the excess free 
energy is a volume integral over the excess free energy 
density, F ex — J dr{§HS + &p} where $jjs is the ex- 
cess free energy density of the reference mixture of hard 
spheres and <i>p is a perturbation contribution due to the 
chain connectivity. 3>hs is evaluated using the elegant 
and inspiring White Bear version 38,39 of Rosenfeld's fun- 
damental measure theory (FMT)iS. Likewise, FMT-style 
weighted densities are used for the polymer contribution, 
$p. We refer the reader to earlier papers for explicit 
formulae^SiSiSi. The equilibrium density profiles were 
found from the condition 



SSI [pp (R) , pc (r )] _ Sn [ PP (R) , pc (r)] 



5p P (R) 
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FIG. 3: The excess polymer segment adsorption isotherm 
Tps for a colloid-polymer mixture with d — 10 and M = 100, 
calculated along a path of constant polymer segment pack- 
ing fraction xps =0.944 727 94. The dashed line marks the 
prewetting transition, while the dotted line denotes bulk co- 
existence. The inset illustrates the nature of the divergence 
of the excess adsorption (open circles) upon approaching bulk 
coexistence. Tps increases linearly as a function of the loga- 
rithm of the undersaturation Apr- The straight solid line is 
a guide to the eye. 



III. RESULTS 

Paricaud et al 30 shown that within the Wertheim ap- 
proach a mixture of colloids and. excluded volume poly- 
mers undergoes a demixing transition into the colloid-rich 
(polymer-poor) and colloid-poor (polymer-rich) phases. 
The properties of the fluid-fluid interface resulting from 
this theory were reported in Refi^&. In the present pa- 
per we investigate the wettability of hard structureless 
walls by such mixtures. The calculations were carried 
out for constant size ratio d = <tc/&ps = 10 and for 
chain lengths M = 80, 100, 120 and 200. This cor- 
responds to q = 0.894, 1, 1.095 and 1.414, and in the 
above we assumed that R g = apM u ops with ap ~ 0.5 
and v — o.5^2iSiiii^. Thus we consider both regimes 
q < 1 and q > 1 . In Fig. ^ we recall the bulk phase 
diagram for the system with M = 100 plotted in the 
polymer segment mole fraction, xps — Pps/iPps + Pc)i 
- total density, pr — pp% + p% , plane. This represen- 
tation proved to be useful in carrying out numerical cal- 
culations presented below. If we disregarded the chain 
connectivity xpg would correspond to the mole fraction 
of the small spheres in a mixture of big and small spheres, 
while ppg and p£) would correspond to the bulk number 
densities of small and big hard spheres, respectively. The 



bulk densities p PS and p c serve as an input to the DFT 
calculations. 

We note that the critical polymer segment mole frac- 
tion for the system depicted in Fig. ^ is 0.976 763 34, 
thus subcritical state points with lower mole fractions 
correspond to the polymer-poor (colloid-rich) side of the 
phase diagram while state points with xps higher than 
this value correspond to the the polymer-rich (colloid- 
poor) side. The numerical procedure of the evaluation of 
the surface phase diagrams relied on monitoring excess 
adsorption isotherms, T a = J dzp a (z) — p& \ a =C, PS, 
and grand potentials (Eq. calculated along paths of 
constant total density (an example of such path labelled 
as "I" is shown in Fig. [I}, 01 along paths of constant poly- 
mer segment mole fraction (the paths marked as "II" and 
"III" in Fig. P. 

The inset to Fig. ^ shows a comparison of the bin- 
odal evaluated for M — 100 resulting from the Werthiem 
TPT1 description with the binodal resulting from the 
free- volume theory for the AOV models for q = 1 . The 
diagrams are plotted in the colloid-polymer packing frac- 
tion representation. For such comparison to be possible 
one has to transform the microscopic polymer segment 
density ppg into the polymer packing fraction defined as 

t)P = 4/3nR 3 g p { pl/M. When compared with the AOV 
model the Wertheim TPT1 approach leads to a shift of 
the colloid and polymer critical packing fractions similar 
to the perturbative treatment of Warren et alr& and the 
geometrical-based approach^. We note that the polymer 
packing fractions are too low when compared with exper- 
imental values^. This is in part due to the fact that, as 
mentioned earlier, the Wertheim TPT1 theory leads to 
the incorrect Flory exponent (y = 0.5 rather than 0.588). 

We start the presentation of the results of DFT cal- 
culations by analyzing the colloid-poor side of the bulk 
phase diagram. The upper panel of Fig. |21 shows the 
excess polymer segment adsorption isotherm evaluated 
along path I, i.e. along a path of constant total density. 
We find the adsorption isotherm to be smooth and small 
in magnitude. We have carefully inspected adsorption 
isotherms calculated along paths similar to the path I 
for the system with M = 100 and other chain lengths 
and have always found similar behavior. Therefore we 
conclude that within the TPT1 formalism mixtures with 
the state points located on the colloid-poor side of the 
phase diagram always stays in the partial wetting/drying 
regime. This is the opposite of what is found in the wet- 
ting studies of the AOV model 26 i 27 i 28 . 

We turn now to the wetting properties of the state 
points located on the polymer-poor side of the phase dia- 
gram. The lower panel of Fig.|21shows the excess polymer 
segment adsorption isotherm calculated along path II, i.e. 
along a path of constant polymer segment mole fraction 
x PS =0.928 889 85. Similar to the path I we find the ad- 
sorption isotherm to be smooth, small in magnitude, and 
finite at bulk coexistence. However, when paths located 
closer to the bulk critical point are considered, a dras- 
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FIG. 4: Average polymer segment (a) and colloid (b) density 
profiles for a colloid-polymer mixture with d = 10 and M = 
100 near a hard wall. The profiles are evaluated for constant 
polymer segment packing fraction xps =0.944 727 94 and for 
total bulk densities pT<?ps = 0.012 3625 (the dashed and solid 
lines correspond to the coexisting thin and thick film phases of 
the prewetting transition, respectively), 0.0126 (dotted lines), 
and 0.012 617 (dash-dotted lines). 



tic change in the behavior of the inhomogeneous colloid- 
polymer mixture is observed. This is illustrated in Fig. [31 
where we show the adsorption isotherm, calculated along 
path III i.e. along a path of constant xps =0.944 727 
94. At lower total bulk densities the excess adsorption 
is very similar to that evaluated along path II. However, 
as pt increases the adsorption isotherm makes a sud- 
den jump (marked by the dashed line) to a large positive 
value and diverges upon approaching bulk coexistence 
(denoted by the vertical dotted line). The nature of this 
divergence is investigated further in the inset to Fig. [3] 
We find that the excess adsorption diverges logarithmi- 
cally and this behavior is characteristic of complete wet- 
ting for short range forces^. We therefore conclude that 
in the TPT1 description the colloid-poor demixed phase 
wets completely the hard wall at coexistence. The wet- 
ting transition is of the first order and is accompanied by 
the prewetting (or thin-thick film) transition. 

Morphology of the demixed mixture close to the hard 
wall is investigated further in Fig. ^ where we show exam- 
ples of the density profiles evaluated at fixed xps =0.944 
727 94 and for several total bulk densities. The dashed 
and solid lines correspond to the coexisting thin and thick 
film profiles (these states are marked as black dots in 
Fig. 01 . The polymer segment profile (cf. Fig.QJi, dashed 
line) corresponding to the thin film phase exhibits oscil- 
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FIG. 5: Surface phase diagram for a colloid polymer mixture 
with d = 10 and M — 100 in the polymer segment mole 
fraction-total density plane. The binodal and the spinodal 
are denoted by the solid and dashed lines, respectively. The 
open circle denotes the bulk critical point. The filled circle 
denotes the wetting point whereas the triangle denotes the 
surface critical point. Thick solid line denotes the prewetting 
line. The inset shows the surface phase diagram in the colloid- 
polymer packing fraction representation. 



lations in the direct proximity of the wall, and the oscil- 
lation period approximately equals the colloid diameter. 
However the excess adsorption is small in magnitude and 
negative. The colloid profile (cf. Fig. 0Jd, dashed line) 
also exhibits oscillations but the structure is much more 
pronounced and the profile exhibits a sharp and high first 
peak with a large contact value. This suggests that in the 
thin film phase the polymers are depleted from the wall, 
whereas the colloids are attracted to the wall and form a 
highly packed first layer. 

On the other hand, in the thick film phase of the poly- 
mer segment density profile (cf. Fig. solid line) one 
observes formation of a thick polymer layer next to the 
wall. As the total density increases (cf. Fig. dotted 
and dash-dotted lines) the polymer-rich layer grows and 
extends over the distance of many colloid diameters. The 
structure of the colloids in the thick film phase is a re- 
verse of the structure of the polymer. The contact value 
of the colloid profile (see the inset to Fig. 0Jd) is reduced 
by more than a half at the prewetting transitional sug- 
gesting that the colloids are depleted from the region ad- 
joining the wall. The colloid-poor region broadens upon 
approaching bulk coexistence (Fig. U>, dotted and dash- 
dotted lines). 

Figures [S] and show the surface phase diagram for 
colloid-polymer mixtures with M = 100 and M = 200, 
respectively. As in simple fluids the prewetting line 
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short-chain fluids^ it was suggested that the wetting 
temperature for a given surface, when divided by the crit- 
ical temperature, is almost independent on M for longer 
chains. Although in the athermal systems considered in 
the present paper the temperature does not have any 
significant role (is a simple scaling parameter), other pa- 
rameters like the reciprocal pressure play the role sim- 
ilar to the temperature in thermal systems. Therefore 
the wetting point pressure rescaled by the critical pres- 
sure P w = P( cr '/P w should barely depend on the chain 
length. The estimated wetting point pressure f3P w a^ for 
systems M = 80, 100, 120, and 200 is 9.75, 7.76, 6.51, 
and 4.25, respectively. The corresponding P w thus are 
0.814, 0.818, 0.824, and 0.835. This agrees nicely with 
the prediction^. 



IV. DISCUSSION 



PS 



FIG. 6: Surface phase diagram for a colloid polymer mixture 
with d = 10 and M — 200 in the polymer segment mole 
fraction-total density plane. The binodal and the spinodal are 
denoted by the solid and dashed lines, respectively. The open 
circle denotes the bulk critical point. The filled circle denotes 
the wetting point whereas the triangle denotes the surface 
critical point. Thick solid line denotes the prewetting line. 
The inset shows the surface phase diagram in the polymer 
segment mole fraction-pressure representation. 



tangentially meets the binodal in the wetting transition 
point (marked by the filled circle) and ends at the sur- 
face critical point (marked by the triangle) . In the poly- 
mer segment mole fraction-total density representation 
the prewetting line is located very close to the binodal 
and the surface critical point is well removed from the 
bulk critical point (marked by the open circle). How- 
ever, when the colloid-polymer packing fraction (cf. the 
inset to Fig. [SJ or the polymer segment mole fraction- 
pressure representation is chosen (cf. the inset to Fig.|BJ), 
the prewetting line is more distant from the binodal. We 
note that for all chain lengths considered the surface crit- 
ical point pressure is lower than the pressure of the bulk 
demixing transition. This curious feature was also found 
in a recent study of the wetting behavior of polymer 
solutions 45 . However in the present paper we do not 
find the lower wetting points and, despite considerable 
effort, we do not find any hints for the layering transi- 
tions in the partial wetting regime. A likely explanation 
to the last observation is that Forsman and Woodward 45 
considered a mixture of polymers and solvent molecules 
with diameter ratio d < 1. Consequently, the individual 
polymer segments were effectively attracted to the wall 
due to the depletion forces induced by the overlap of the 
excluded volumes of the polymer segments and the wall. 
In the present study this effect does not come into play. 
In a recent study of wettability of solid surfaces by 



That the rescaled P w barely depends on M has a sim- 
ple explanation. As demonstrated in Refsi^SiSi, in the 
limit of an infinite chain length M — ► oo, and for fixed 
d, the bulk critical pressure (and also the colloid criti- 
cal packing fraction) tend to a finite, nonzero value with 
the leading order correction ~ 0(l/vM). This is the so- 
called "protein limit" of colloid-polymer mixtures^. But 
then the wetting point pressure also should tend to a con- 
stant in that limit. Therefore the rescaled wetting point 
pressure should barely depend on M for longer chains. 
We also note that the systematic deviations of P w from a 
constant suggest that the leading order chain-length de- 
pendent correction to P w is of different order than that 
for the bulk critical point. However the precise form of 
these corrections is rather difficult to estimate and we 
leave this problem for future work. 

Wetting in the AOV model for q < 1 can be readily 
explained by invoking depletion interactions. The over- 
lap of excluded volumes of a colloid and a wall is larger 
than the overlap between two colloidal particles and this 
gives rise to net colloid-wall attraction. Consequently the 
colloid-rich phase favors the hard wall. Our results are 
in opposition to the AOV model results and show some 
features that are characteristic of wetting in athermal 
polymer solutions. 

The suggestion that the polymer-rich phase should wet 
the hard wall seems reasonable when q > 1. However, we 
find this behavior also for q < 1 i.e. for d = 10, M = 80 
(in order to be completely assured that our calculations 
are representative of the TPT1 description and are not 
a result of e.g. an artifact connected with the incorrect 
scaling of the polymer gyration radius, we checked that 
even for colloid-to-polymer-segment size ratios as large as 
d = 20 and for chains as short as M = 40 the polymer- 
rich phase wets completely the hard wall). This is at vari- 
ance with experimental results. We connect this discrep- 
ancy with the inaccurate description of the bulk phase 
behavior resulting from the TPT1 theory, and identify 
two possible sources of errors. The first issue is con- 
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nected with the incorrect scaling of the gyration radius 
and of the polymer second virial coefficient with M. The 
consequence is that the dilute and semidilute regimes of 
long chains are not described accurately. This problem 
has already been recognized in the literature. Vega et al. 
proposed a simple correction that significantly improves 
the accuracy of the TPT1 description in the low density 
limfbali. A similar treatment, which combines the field- 
theoretic approach with liquid-state theory, was proposed 
by Luc 42 . However, at present it is not clear how these 
ideas can be extended to mixtures. To explain the sec- 
ond deficiency in the TPTf description of colloid-polymer 
mixtures we recall that the Wertheim approach, as any 
other perturbation theory, relies on the good description 
of the reference system. In this case it is a mixture of big 
and small hard spheres and we used the usual Boublik- 
Mansoori-Carnahan-Starling-Leland (BMCSL) equation 
of stat o 47 i 48 . However, it was argued— that for highly 
asymmetric hard sphere mixtures a fluid-fluid phase sep- 
aration occurs, which is metastable to a broad fluid-solid 
phase transition. Consequently the BMCSL equation of 
state is not very accurate. The reference system in the 
TPTf approach to colloid-polymer mixtures consists of a 
very low packing fraction of small spheres, therefore the 
potential impact of this problem seems to be moderate. 

V. CONCLUSIONS 

We have used a microscopic density functional theory 
based on Wertheim's first order thermodynamic pertur- 



* Electronic address: pawel@paco.umcs.lublin.pl 

1 W.C.K. Poon, J. Phys.: Condens. Matter 14, R859 (2002). 

2 R. Tuinier, J. Rieger, and C.G. de Kruif, Adv. Colloid 
Interface Sci. 103, 1 (2003). 

3 A.P. Cast, C.K. Hall, and W.B. Russell, J. Colloid Inter- 
face Sci. 96, 251 (1983). 

4 S. Asakura and F. Oosawa, J. Chem. Phys. 22, 1255 
(1954). 

5 A. Vrij, Pure Appl. Chem. 48, 471 (1976). 

6 H.N.W. Lekkerkerker, W.C.K. Poon, P.N. Pusey, A. 
Stroobants, and P.B. Warren, Europhys. Lett. 20, 559 
(1992). 

7 J.M. Brader, R. Evans, and M. Schmidt, Mol. Phys. 101, 
3349 (2003). 

8 M. Fuchs and K.S. Schweizer, J. Phys.: Condens. Matter 
14, R239 (2002). 

9 M. Fuchs and K.S. Schweizer, Europhys. Lett. 51, 621 
(2000). 

10 S.A. Shah, Y.L. Chen, K.S. Schweizer, and C.F. Zukoski, 
J. Chem. Phys. 118, 3350 (2003). 

11 P.G. Boluis, A. A. Louis, and J.-P. Hansen, Phys. Rev. Lett. 
89, 128302 (2002). 

12 P.B. Warren, S.M. Ilett, and W.C.K. Poon, Phys. Rev. E 
52, 5205 (1995). 

13 D.G.A.L. Aarts, R. Tuinier, and H.N.W. Lekkerkerker, J. 



bation theory to study the wetting behavior of athermal 
mixtures of colloids and excluded- volume polymers. In 
opposition to the wetting behavior in the AOV model we 
have found the polymer-rich phase to wet a hard wall. 
The wetting transition is of first order and is accompa- 
nied by the prewetting transition. We have not found 
any hints for the layering transitions in the partial wet- 
ting regime. The rescaled wetting point pressure is very 
similar for mixtures with different chain lengths and the 
surface critical point pressure is lower than the pressure 
of the bulk demixing transition. We have found our re- 
sults to resemble the wetting behavior in athermal poly- 
mer solutions. 

Our considerations show that it is certainly a challenge 
to devise a monomer-resolved theory for colloid-polymer 
mixtures that is accurate in both the 'colloid', and 'pro- 
tein' limits. However, we feel that the Wertheim ap- 
proach could be a good candidate but this theory should 
incorporate the correct scaling behavior in the dilute 
polymer regime and an accurate description of the ref- 
erence system. 



Acknowledgments 



This work has been supported by KBN of Poland under 
the Grant 1P03B03326 (years 2004-2006). 



Phys.: Condens. Matter 14, 7551 (2002). 

14 M. Schmidt, A.R. Denton, and J.M. Brader, J. Chem. 
Phys. 118, 1541 (2003). 

15 A. Moncho-Jorda, J. Dziubiella, J.P.Hansen, and A. A. 
Louis, J. Phys. Chem. B 109, 6640 (2005). 

16 P. Bolhuis, E.J. Meijer, and A.A. Louis, Phys. Rev. Lett. 
90, 68304 (2003). 

17 A. Moncho-Jorda, A.A. Louis, P.G. Bolhuis, and R. Roth, 
J. Phys.: Condens. Matter 15, S3429 (2003). 

18 D.G.A.L. Aarts, M. Schmidt, and H.N.W. Lekkerkerker, 
Science 304, 847 (2004). 

19 D.G.A.L. Aarts, J. Phys. Chem. B 109, 7407 (2005). 

20 E.H.A. de Hoog, H.N.W. Lekkerkerker, J. Schulz, and G.H. 
Findenegg, J. Phys. Chem. B 103, 10657 (1999). 

21 D.G.A.L. Aarts, J.H. van der Wiel, and H.N.W. Lekkerk- 
erker, J. Phys.: Condens. Matter 15, S245 (2003). 

22 B.-H. Chen, B. Payandeh, and M. Robert, Phys. Rev. E 
62, 2369 (2000). 

23 B.-H. Chen, B. Payandeh, and M. Robert, Phys. Rev. E 
64, 042401 (2001). 

24 W.K. Wijting, N.A.M. Besseling, and M.A. Cohen Stuart, 
Phys. Rev. Lett. 90, 196101 (2003). 

25 W.K. Wijting, N.A.M. Besseling, and M.A. Cohen Stuart, 
J. Phys. Chem. B 107, 10565 (2003). 

26 J.M. Brader, R. Evans, M. Schmidt, and H. Lowen, J. 



7 



Phys.: Condens. Matter 14, LI (2002). 

27 M. Dijkstra and R. van Roij, Phys. Rev. Lett. 89, 208303 
(2002). 

28 D.G.A.L. Aarts, R.P.A. Dullens, H.N.W. Lekkerkerker, D. 
Bonn, and R. van Roij, J. Chem. Phys. 120, 1973 (2004). 

29 D.G.A.L. Aarts and H.N.W. Lekkerkerker, J. Phys.: Con- 
dens. Matter 16, S4231 (2004). 

30 P. Paricaud, S. Varga, and G. Jackson, J. Chem. Phys. 
118, 8525 (2003). 

31 M.S. Wertheim, J. Chem. Phys. 87, 7323 (1987). 

32 C. Vega and L.G. MacDowell, J. Chem. Phys. 114, 10411 
(2001). 

33 T. Boublik, Mol. Phys. 68, 191 (1989). 

34 T. Boublik, C. Vega, and M. Diaz Pena, J. Chem. Phys. 
93, 730 (1990). 

35 C. Vega, J.M. Labaig, L.G. MacDowell, and E. Sanz, J. 
Chem. Phys. 113, 10398 (2000). 

36 P.Bryk, J. Chem. Phys. 122, 064902 (2005). 

37 Y.-X. Yu and J. Wu, J. Chem. Phys. 117, 2368 (2002). 

38 R. Roth, R. Evans, A. Lang, and G. Kahl, J. Phys.: Con- 
dens. Matter 14, 12063 (2002). 

39 Y.-X. Yu and J. Wu, J. Chem. Phys. 117, 10156 (2002). 



40 Y. Rosenfeld, Phys. Rev. Lett. 63, 980 (1989). 

41 J. Dautenhahn and C.K. Hall, Macromolecules 27, 5399 
(1994). 

42 L. Lue, J. Chem. Phys. 112, 3422 (2000). 

43 S. Dietrich, in Phase Transitions and Critical Phenomena, 
ed. C. Domb and J.L. Lebowitz, Academic, London, 1988, 
Vol. 12, p. 1. 

44 The sum of the contact values of the colloid and the av- 
erage polymer segment density profiles is the same for the 
coexisting thin and thick film phases as dictated by the 
appropriate sum rule. 

45 J. Forsman and C.E Woodward, Phys. Rev. Lett. 94, 
118301 (2005). 

46 P. Bryk and S. Sokolowski, J. Chem. Phys. 121, 11314 
(2004). 

47 T. Boublik, J. Chem. Phys. 53, 471 (1970). 

48 G.A. Mansoori, N.F. Carnahan, K.E. Starling, and T.W. 
Leland, J. Chem. Phys. 54, 1523 (1971). 

49 M. Dijkstra, R. van Roij, and R. Evans, Phys. Rev. E 59, 
5744 (1999). 



